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ABSTRACT 



The transiting extrasolar planet XO-3b is remarkable, with a high mass and eccentric orbit. These unusual characteristics make it 
interesting to test whether its orbital plane is parallel to the equator of its host star, as it is observed for other transiting planets. 
We performed radial velocity measurements of XO-3 with the SOPHIE spectrograph at the 1.93-m telescope of Haute-Provence 
Observatory during a planetary transit, and at other orbital phases. This allowed us to observe the Rossiter-McLaughlin effect and, 
together with a new analysis of the transit light curve, to refine the parameters of the planet. The unusual shape of the radial velocity 
anomaly during the transit provides a hint for a nearly transverse Rossiter-McLaughlin effect. The sky-projected angle between the 
planetary orbital axis and the stellar rotation axis should be A = 70° ± 15° to be compatible with our observations. This suggests that 
some close-in planets might result from gravitational interaction between planets and/or stars rather than migration due to interaction 
with the accretion disk. This surprising result requires confirmation by additional observations, especially at lower airmass, to fully 
exclude the possibility that the signal is due to systematic effects. 

Key words. Techniques: radial velocities - Stars: individual: GSC03727-01064 - Stars: planetary systems: individual: XO-3b 



1. Introduction 



Johns-Krull et al. (2008 ) announced the detection of XO-3b, an 
extra-solar planet transiting its F5V parent star with a 3.2-day or- 
bital period. Transiting planets are of particular interest as they 
allow measurements of parameters including orbital inclination 
and planet radius, mass and density. Moreover, follow-up ob- 
servations can also be performed during transits or anti-transits, 
yielding physical constraints on planetary atmospheres. 

Among the forty transiting extra-solar planets known to date, 
XO-3b is particular as it is among the few on an eccentric or- 
bit, together with HP 147506b (Bakoset al. [20071), HD 17156b 
(Fisc her et al. [20071 Bar bieri et al. p2007l and GJ436b (Butler 
et al. 120041 Gillon et al. [2QQ7l >. XO-3b is also the most mas- 
sive transiting planet known to date. Most of the sixty known 
extrasolar planets, with and without transits, with orbital peri- 
ods shorter than five days have masses below 2M Jup ; XO-3b 
is actually one of the rare massive close-in planets. It is just 
at the limit between low-mass brown dwarfs and massive plan- 
ets, 13M Jup , which is defined by the deuterium burning limit. 
There was a quite large uncertainty on the planetary parame- 
ters of XO-3b and its host star. Indeed, Johns-Krull et al. (2008 ) 
presented a spectroscopic analysis favoring large masses and 
radii (M p 13.25 M Jup , R p 1.95 R Jup , M* 1.41 M G , and 
R+ - 2.13 R ), whereas their light curve analysis suggests lower 



* Based on observations collected with the SOPHIE spectrograph 
on the 1.93-m telescope at Observatoire de Haute-Provence (CNRS), 
France, by the SOPHIE Consortium (program 07A.PNP.CONS). 



values (M p 12.03 M Jup , R p 1.25 R Jup , M* 1. 24 M Q and 
s 1.48 Rq) [see however Sect. [5] and Winn et al. (I2008ab l. 
The fast rotati ng st ar XO-3 (Vsin/ = 18.5 kms" 1 ; 
Johns-Krull et al. 2008) is a favorable object for Rossiter- 
McLaughlin effect observations. This effect (Rossiter 1924; 
McLaughlin 1 1 924t occurs when an object transits in front of a 
rotating star, causing a distortion of the stellar lines profile, and 
thus an apparent anomaly in the measured radial velocity of the 
star. The shape of the disturbed radial velocity curve allows one 
to determine whether the planet is orbiting in the same direction 
as its host star is rotating, and more generally to measure the 
sky-projected angle between the planetary orbital axis and the 
stellar rotation axis, usually noted A (see, e.g., Ohta et al. 2005; 
Gimenez 2006a; Gaudi & Winn 2007 ). A stellar spin axis not 
aligned with the orbital angular momentum of a planet (A + 0°) 
could reflect processes in the planet formation and migration, 
or interactions with perturbing bodies (see, e.g., Malmberg et 
al. 120071 Chatterjee et al. 120071 Nagasawa et al. 120081 ). Solar 
System asteroids are examples of objects whose orbital axes can 
be misaligned from the Sun spin axis by over 30°. 

Up to now, spectroscopic transits have been detected 
for eight exoplanets: HD 209458b (Queloz et al. 120001 , 
HD 189733b (Winn et al.l2006l HD 149026b (Wolf et al. 120071, 
TrES-1 (Narit a et a l. 120071 ). HD 147506b (Wi nn et al. [20071 
Loeillet et al. [JOOHJt, HD 171 56b (Narita et al. 120081 CoRoT - 
Exo-2b (Bouchy et al. 120081 ) and TrES-2 (Winn et al. H)08b ). 
For all of these targets the stellar rotation is prograde relative 
to the planet orbit, and the sky-projected A angle is close to 
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zero for most of them. So the axes of the stellar spins are prob- 
ably parallel to the orbital axes, as expected for planets that 
formed in a protoplanetary disc far from the star and that later 
migrated closer-in. Three systems have error bars on the A angle 
that do not include 0°: TrES-1 (A = 30° ±21°), CoRoT-Exo- 
2b (A = 7.2 ± 4.5°) and HD 17156b (A = 62° ± 25°). However, 
those cases have the largest error bars on A, and no firm detection 
of misalignment has yet been claimed. Barbieri et al. (2008]' re- 
cently presented new radial velocity measurements of HD 17156 
secured during a transit, which agree with a spin-orbit alignment. 

Approximate spin-orbit alignment therefore seems typical 
for exoplanets, as it is for planets in the Solar System. The un- 
usual parameters of XO-3b make a test of whether it agrees with 
this apparent behavior interesting. We present here new mea- 
surements of XO-3 radial velocity performed during a transit 
and at other orbital phases. These data refine the orbital param- 
eters and provide a hint of detection for a transverse Rossiter- 
McLaughlin effect, i.e. a A angle possibly near 90°. We also 
present a revised analysis of the transit light curve. 

2. Observations 

We observed the host star XO-3 (GSC 03727-01064, m v = 9.91) 
with the SOPHIE instrument at the 1 .93-m telescope of Haute- 
Provence Observatory, France. SOPHIE is a cross-dispersed, en- 
vironmentally stabilized echelle spectrograph dedicated to high- 
precision radial velocity measurements (Bouchy et al. 2006). We 
used the high-resolution mode (resolution power R = 75, 000) 
of the spectrograph, and the fast-read-out-time mode of the 
4096 x 2048 15-yum-pixel CCD detector. The two optical-fiber 
circular apertures were used; the first one was centered on the 
target, and the second one was on the sky to simultaneously mea- 
sure its background. This second aperture, 2' away from the first 
one, was used to estimate the spectral pollution due to the moon- 
light, which can be quite significant in these 3"-wide apertures 
(see Sect. [3]l. 

We acquired 36 spectra of XO-3 during the night of January 
28th, 2008 (barycentric Julian date BJD = 2454494.5), 
where a full coverage of the planetary transit was observed. 
Another 19 spectra were acquired at other orbital phases dur- 
ing the following two months. Table [T] summaries the 55 spectra 
finally acquired. 

The exposure times range from 6 to 30 minutes in order to 
reach as constant the signal-to-noise ratio as possible. Indeed, 
SOPHIE radial velocity measurements are currently affected by 
a systematic effect at low signal-to-noise ratio, which is proba- 
bly due to CCD charge transfer inefficiency that increases at low 
flux level. A constant signal-to-noise ratio through a sequence 
of observations reduces this uncertainty. The different exposure 
times needed to reach similar signal-to-noise ratios reflect the 
variable throughputs obtained, due to various atmospheric con- 
ditions (seeing, thin clouds, atmospheric dispersion). The sky 
was clear on the night whom the XO-3b transit was observed, 
but the airmass ranged from 1.2 to 3.1 during this ~6-hour ob- 
servation sequence; the exposure times therefore increased dur- 
ing the transit observation. They remain short enough to provide 
a good time sampling (20 measurements during the ~3 hours of 
the transit). The 19 measurements outside the transit night were 
performed at airmasses better than 1 .4 but with conditions vary- 
ing from photometric to cloudy. 

Exposures of a thorium-argon lamp were performed every 
2-3 hours during each observing night. Over 2-3 hours, the ob- 
served drifts were typically ~ 3 m s -1 , which is thus the accuracy 
of the wavelength calibration of our XO-3 SOPHIE spectra; this 



Table 1. Radial velocities of XO-3 measured with SOPHIE. 



BJD 


RV 


±lcr 


exp. time 


S/N p. pix. 


-2400 000 


(km s" 1 ) 


(kms- 1 ) 


(sec) 


(at 550 nm) 


Planetary transit 










54494.4461 


-11.298 


0.020 


600 


54 


54494.4526 


-11.300 


0.026 


403 


42 


54494.4578 


-11.345 


0.027 


373 


40 


54494.4625 


-11.379 


0.027 


370 


40 


54494.4675 


-11.390 


0.028 


370 


39 


54494.4721 


-11.354 


0.028 


370 


39 


54494.4767 


-11.460 


0.029 


370 


38 


54494.4813 


-11.448 


0.028 


381 


39 


54494.4861 


-11.419 


0.027 


380 


40 


54494.4913 


-11.411 


0.027 


380 


40 


54494.4960 


-11.457 


0.028 


380 


39 


54494.5007 


-11.463 


0.028 


380 


39 


54494.5054 


-11.549 


0.028 


380 


39 


54494.5101 


-11.482 


0.028 


380 


40 


54494.5152 


-11.589 


0.028 


380 


39 


54494.5202 


-11.652 


0.028 


434 


39 


54494.5254 


-11.611 


0.027 


403 


40 


54494.5303 


-11.679 


0.027 


392 


40 


54494.5352 


-11.691 


0.036 


395 


40 


54494.5413 


-11.667 


0.035 


509 


43 


54494.5474 


-11.782 


0.034 


500 


44 


54494.5535 


-11.700 


0.034 


500 


42 


54494.5599 


-11.649 


0.035 


554 


44 


54494.5665 


-11.785 


0.034 


531 


44 


54494.5738 


-11.811 


0.036 


530 


44 


54494.5806 


-11.783 


0.033 


602 


45 


54494.5880 


-11.867 


0.033 


622 


44 


54494.5957 


-11.767 


0.033 


653 


45 


54494.6036 


-11.785 


0.035 


650 


46 


54494.6123 


-11.740 


0.034 


682 


45 


54494.6210 


-11.783 


0.035 


775 


45 


54494.6305 


-11.880 


0.033 


801 


45 


54494.6407 


-11.843 


0.034 


906 


44 


54494.6532 


-11.920 


0.036 


964 


42 


54494.6668 


-11.964 


0.037 


1321 


43 


54494.6822 


-11.913 


0.049 


1300 


38 


Other orbital phases: 








54496.2649 


-12.723 


0.050 


1202 


22 


54497.2609 


-10.156 


0.029 


655 


37 


54499.2765 


-13.006 


0.030 


1003 


36 


54501.2926 


-12.433 


0.031 


775 


34 


54501.4628 


-12.756 


0.033 


775 


34 


54502.2730 


-13.068 


0.024 


614 


44 


54503.2614 


-10.936 


0.030 


907 


35 


54503.4700 


-10.182 


0.040 


1806 


29 


54504.4321 


-12.398 


0.028 


645 


39 


54505.2889 


-13.132 


0.023 


635 


46 


54506.2904 


-11.593 


0.025 


755 


43 


54511.4534 


-13.041 


0.038 


600 


48 


54512.4618 


-12.246 


0.036 


600 


41 


54513.3091 


-10.360 


0.052 


600 


48 


54516.3517 


-10.176 


0.046 


600 


39 


54516.4540 


-10.267 


0.071 


802 


32 


54551.3044 


-10.316 


0.032 


1274 


37 


54553.3002 


-13.135 


0.033 


999 


36 


54554.3114 


-11.004 


0.019 


999 


57 



is good enough for the expected signal. We did not use simulta- 
neous calibration to keep the second aperture available for sky 
background estimation. During the night of January 28th, 2008, 
we performed a thorium-argon exposure before the transit and 
another one after the sequence, about six hours later. The mea- 
sured drift was particularly low this night, 1 ms _1 in six hours, 
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which makes us confident the wavelength calibration did not un- 
expectedly drift during the observation of the XO-3b transit. 



3. Data reduction 

We extracted the spectra from the detector images and mea- 
sured the radial velocities using the SOPHIE pipeline. Following 
the techniques described by Baranne et al. (1996) and Pepe et 
al. J2002), the radial velocities were obtained from a weighted 
cross-correlation of the spectra with a numerical mask. We used 
a standard G2 mask constructed from the Sun spectrum atlas 
including more than 3500 lines, which is well adapted to the 
F5V star XO-3. We eliminated the first eight spectral orders of 
the 39 available ones from the cross-correlation; these blue or- 
ders are particularly noisy, especially for the spectra obtained at 
the end of the transit, when the airmass was high. The resulting 
cross-correlation functions (CCFs) were fitted by Gaussians to 
get the radial velocities, as well as the width of the CCFs and 
their contrast with respect to the continuum. The uncertainty on 
the radial velocity was computed from the width and contrast 
of the CCF and the signal-to-noise ratio, using the empirical re- 
lation detailed by Bouchy et al. ([2005 ) and Collier Cameron et 
al. (12007b . It was typically around 25 ms _I during the night of 
the transit, and between 20 and 45 ms~' the remaining nights. 
The large Vsin/ of this rotating stars makes the uncertainty 
slightly larger than what is usually obtained for such signal-to- 
noise ratios with SOPHIE. 

Some measurements were contaminated by the sky back- 
ground, including mainly the moonlight. As the G2 mask 
matches the XO-3 spectrum as well as the Sun spectrum re- 
flected by the Moon and the Earth atmosphere, the moonlight 
contamination can distort the shape of the CCF and thus shift the 
measured radial velocity. During our observations, the 29-km s _I 
wide (FWHM) CCF of XO-3 is at radial velocities between -13 
and -10 kms -1 , whereas the moonlight was centered near the 
barycentric Earth radial velocity, between -23 and -20 kms -1 . 
Thus moonlight contamination tends to blueshift the measured 
radial velocities. Following the method described in Pollacco et 
al. ([2008 ) and Barge et al. ((2008 ), we estimated the Moon con- 
tamination thanks to the second aperture, targeted on the sky, and 
then subtracted the sky CCF from the star CCF (after scaling by 
the throughput of the two fibers). Five exposures with too strong 
contamination were not used. We estimated the accuracy of this 
correction on one hand by correcting uncontaminated spectra, 
and on the other hand by correcting uncontaminated spectra on 
which we have added moonlight contaminations. Comparisons 
of the corrected velocities to the uncontaminated ones show that 
the method works well up to ~ 500 ms~' shifts, with an uncer- 
tainty of 1/9 of the correction to which a minimum uncertainty 
of 25 ms _1 is quadratically added. 

The second half of the transit night measurements was con- 
taminated by moonlight, with sky CCFs contrasted between 2 
and 5 % of the continuum, whereas the XO-3 CCF has a contrast 
of 8 %. This implied sky corrections < 150 m s _1 (except for the 
very last exposure where it was ~ 300 ms _1 ), with uncertainties 
in the range 25-30 m s _1 (40 m s _1 for the last exposure). Five ex- 
posures obtained later at different phases were contaminated by 
the moonlight; corrections of 100 to 500 ms~' were computed, 
with uncertainties in the range 30-60 m s . 

The final radial velocities are given in Table[T]and displayed 
in Figs. [T] and [2] The error bars are the quadratic sums of the 
different error sources (photon noise, wavelength calibration and 
drift, moonlight correction). 



XO-3 SOPHIE 




95 100 105 110 115 120 125 130 135 140 145 150 155 
BJD-2454400 [days] 



Fig. 1. Top: Radial velocity measurements of XO-3 as a func- 
tion of time, and Keplerian fit to the data (without transit). Only 
the 23 measurements used for the fit are displayed. The or- 
bital parameters corresponding to this fit are reported in Table [2] 
Bottom: Residuals of the fit with 1-cr error bars. 



4. Determination of the planetary system 
parameters 

4.1. Refined orbit 

The radial velocities measurements presented by Johns-Krull et 
al. ([2008 ) have a typical accuracy of ~ 160 m s _1 . Those secured 
with SOPHIE are about five times more accurate, so they allow 
for a refinement of the original parameters of the system. We 
made a Keplerian fit of the first four SOPHIE measurements per- 
formed during the transit night and those performed afterwards, 
at other phases, first using the orbital period from Johns-Krull et 
al. (12008). For the refinement of the orbit we did not use most 
of the data secured during the night of January 28th, 2008, in 
order to remain free from alteration due to transit anomalies and 
possible systematic effects due to large-airmass observations. 

The standard deviation of the residuals to the fit is cr(0-C) = 
29 ms _I , implying a^ 2 of 15.3, which is acceptable according 
the low degrees of freedom, v = 18. The 29 ms _1 dispersion of 
the measurements around the fit is similar to the errors on the in- 
dividual radial velocity measurements; these estimated error bars 
thus are approximatively correct. The fits are plotted in Figs. [T] 
and[2j the derived orbital parameters are reported in Table [2] to- 
gether with error bars, which were computed from^ 2 variations 
and Monte Carlo experiments. They agree with the Johns-Krull 
et al. (120081 ) parameters but the error bars are reduced by fac- 
tors of three to six. The largest difference is on the eccentricity, 
which we found 1.6 cr larger than Johns-Krull et al. (2008 ). The 
residuals, plotted as a function of time in the bottom panel of 
Fig.[T] do not show any trend that might suggest the presence of 
another companion in the system over two months. 
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Fig. 2. Phase-folded radial velocity measurements of XO-3 (cor- 
rected from the velocity V r = -12.045 km s _1 ) as a function of 
the orbital phase, and Keplerian fit to the data. Orbital parame- 
ters corresponding to this fit are reported in Table [2] For display 
purpose, all the measurements performed during the transit night 
are plotted here. However, only the first four measurements of 
the transit night are used for the orbit fit, together with 19 mea- 
surements secured at other orbital phases (see § 4.1 1. Figs.[5]and 
[TJdisplay a magnification on the transit night measurements. 



Table 2. Fitted orbit and planetary parameters for XO-3b. 



Parameters 


Values and 1-cr error bars 


Unit 


V r 


-12.045 + 0.006 


km s -1 


P 


3.19161 ±0.00014 


days 


e 


0.287 ± 0.005 






-11.3 ± 1.5 




K 


1.503 ±0.010 


kms -1 


To (periastron) 


2454493.944 ±0.009 


BJD 


cr(0 - C) 


29 


ms~' 


reduced x 1 


0.85 




N 


23 




t c (transit) 


2 454 494.549 ±0.014 


BJD 


M* 


1.3 ±0.2 


M Q 


R* 


1.6 ±0.2 


Ro 


Mp sin i 


12.4 ± 1.9 f 


M Jup 


i 


82.5 ± 1.5 


o 


M p 


12.5 ± 1.9 f 


M Jup 


R P 


1.5 ±0.2 


Rjup 


A 


70 ± 15 





and is used for the fits plotted in Figs.[T]and[2] Adding HJS and 
HET data does not significantly change the other orbital parame- 
ters nor their uncertainties. For the global fit using the radial ve- 
locities from the three instruments, we did not use the last HET 
measurement, performed during a transit (see Sect. [6J. 

The Keplerian fit of the new SOPHIE radial velocity mea- 
surements also improves the transit ephemeris, as the photomet- 
ric transits reported by Johns-Rrull et al. (2008 ) were secured 
between December 2003 and March 2007, one hundred or more 
XO-3b revolutions before the January 28th, 2008 transit. The 
mid-point of this transit predicted from the Keplerian fit of the 
SOPHIE radial velocity measurements is t c = 2454494.549 + 
0.014 (BJD), i.e. just a few minutes earlier than the prediction 
from Johns-Rrull et al. (2008 ). The uncertainty on this transit 
mid-point is +20 minutes (or ±0.004 in orbital phase). 

In order to reduce this uncertainty, we observed a recent 
photometric transit of XO-3b with a 30-cm telescope at the 
Teide Observatory, Tenerife, Spain, on February 29th, 2008 
(Fig. [3 1. Weather conditions were poor and we therefore ana- 
lyzed the transit with a fixed model based on the algorithm of 
Gimenez (2006b ). The fixed parameters were the ratio between 
the radii of the star and of the planet k = 0.0928, the sum of the 
projected radii rr = 0.2275, the inclination i = 79.3°, and the ec- 
centricity e = 0.26. We then scanned different mid-transit times 
and found t c = 2454526.4668 + 0.0026 (BJD) from varia- 
tions. This reflects photon noise only; fluctuations due to poor 
weather may introduce additional uncertainties. By taking into 
account for the uncertainty on the orbital period, this translates 
into t c = 2454494.5507 ± 0.0030 (BJD) for the spectroscopic 
transit that we observed with SOPHIE on January 28th, 2008, i.e. 
ten revolutions earlier. That is just two minutes after the above 
prediction from SOPHIE ephemeris, and the uncertainty on this 
transit mid-point is +4.3 minutes (or ±0.0009 in orbital phase). 



-0.01 



o.oo - 



0.01 



f: using M* = 1.3 ± 0.2 M G 



Fitted alone, the 23 SOPHIE measurements have too short 
time span (60 days) to measure the period more accurately 
than Johns-Rrull et al. (2008) from photometric observations 
of twenty transits. A 1.5-year time span is obtained when the 
SOPHIE measurements are fitted together with the radial ve- 
locities measured by Johns-Rrull et al. (2008 ) using the tele- 
scopes Harlan J. Smith (HJS) and Hobby-Eberly (HET). This 
longer time span allows a more accurate period measurement. 
We obtained P = 3.19168 ± 0.00015 days from the fit using 
the three datasets, in agreement with the photometric one, and 
with a similar uncertainty. The final period reported in Table [2] 
(P = 3.19161 ± 0.00014 days) reflect these two measurements 



0.02 L_ 
0.30 



0.35 



0.40 



0.45 0.50 

HJD - 2434526 



0.55 0.60 0.65 



Fig. 3. Light curve of XO-3 observed at the Teide Observatory, 
Tenerife, during the transit of February 29th, 2008. The tran- 
sit fit (solid line) provides t c = 2454526.4668 + 0.0026 = 
2454494.5507 ± 0.0030 (BJD). 
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4.2. Transit light curve fit revisited 



Johns-Krull et al. (2008 ) point out that the host star radius ob- 
tained from the spectroscopic parameters (temperature, grav- 
ity, metallicity) combined with stellar evolution models, 
2.13 R Q , is incompatible with the value obtained from the shape 
of the transit light curve, namely R+ 1.48 R Q . Indeed, a large 
stellar radius implies a large planetary radius (to account for the 
depth of the transit) and a large inclination angle (to account 
for the duration of the transit), but the time from the first to the 
second contacts (ingress) and third to fourth contacts (egress) 
predicted for such an inclination are too long when compared 
to the observed transit light curve (see the upper panel of Fig. 9 
in Johns-Krull et al. 120081 ). Formal uncertainties on the stellar 
spectroscopic parameters and the photometric measurements are 
insufficient to account for the mismatch. Since there can be only 
one value of the real stellar radius, this must be due to systematic 
uncertainties on the spectroscopic parameters, or the parameter 
derivation from the photometric data, or both. We revisit these 
analyses below, using the photometric data from Johns-Krull et 
al. (f2008 ) and the parameters of the Keplerian orbit obtained in 
§ 4. 1 from the SOPHIE radial velocity measurements. 



Regarding the spectroscopic parameters, the formal uncer- 
tainties stated by Johns-Krull et al. ( 120081 ), e.g. 0.06 dex for the 
gravity log g or 0.03 dex for the metallicity Z, are particularly 
small. Since these are used in combination with stellar evolu- 
tion models, even if the actual uncertainties on the observations 
are small, systematic uncertainties are known to be present in 
the models themselves. Also, precise gravity measurements are 
difficult to obtain from stellar spectra. We therefore set a floor 
level of effective uncertainties in the confrontation with stellar 
evolution models of 100 K in temperature, of 0. 1 dex in logg, or 
0.1 dex in Z (see e.g. discussion in Santos et al. 2004 and Pont 
& Ever l2004l ). 

Regarding the photometric data, we estimated the uncertain- 
ties including systematics effects with "segmented bootstrap" 
analysis (Jenkins et al. [2002 ; Moutou et al. 120041) . According 
to Pont et al. (2006), correlated noise usually dominates the to- 
tal parameter error budget for ground-based transit light curves. 
The segmented bootstrap consist of repeating the fit on realiza- 
tions of the data with individual nights selected at random. The 
photometric follow-up for XO-3 by Johns-Krull et al. (2008) 
consists of ten individual nights. Since the sequencing of the 
data within each night is preserved, this method provide error 
estimates that takes into account the actual correlated noise in 
the data. We find much larger uncertainties on the impact pa- 
rameter than the photon-noise uncertainties. This is corrobo- 
rated by the discussion in Bakos et al. (2006 ) of the case of 
HD 189733. With a much deeper transits and a similar number of 
high-precision photometry transits covered from several obser- 
vatories, they found that the determination of the stellar radius 
from the photometric data produced an error of ~ 15 %, consis- 
tent with the discussion in Pont et al. (2006 ). 

To estimate the probability distribution of the radius of 
XO-3 given the available photometric and spectroscopic ob- 
servations, and the a priori assumption that the star is located 
near theoretical stellar evolution tracks, we use a Bayesian ap- 
proach. As discussed in Pont & Eyer (2004), such approach is 
needed for realistic parameter estimates when the uncertainties 
are not small compared to the total parameter space and the re- 
lation between parameters and observable quantities are highly 
non-linear. We thus calculate the posterior probability distribu- 
tion of the stellar radius R+, using Bayes' theorem with stel- 
lar evolution models and prior probability distributions suitable 
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Fig. 4. Posterior probability distribution function for the stellar 
radius of XO-3 obtained from Bayesian approach. Dashed line: 
Using only the constraints from the light curves of Johns-Krull 



et al. ( 120081 ) and the parameters of the Keplerian orbit (§ 4. 1 



Dotted line: Using only the constraints from the spectroscopic 
parameters. Solid line: Using all these constraints together. 



for a Solar-Neighbourhood magnitude-limited sample, as dis- 
cussed in Pont & Eyer (2004) in the context of the Geneva- 
Copenhaguen survey. 

The posterior probability distribution for the stellar radius is 
calculated according to Bayes' theorem: 



P(R*\PS) 



S 



where R+ is the stellar radius, P the photometric observations 
and S the spectroscopic observations. The first two terms on 
the right are the likelihood of the photometric and spectroscopic 
observations, exp(-l/2;^ 2 ), the last term is the a priori distri- 
bution of R+. The integral covers the mass, age and metallic- 
ity parameters. The stellar evolution models provide the func- 
tion R+ = R+(M*,Z, age). For more detailed explanations of the 
method see Pont & Eyer J2004I) . 

Fig.|4]displays the posterior probability distribution function 
for the stellar radius obtained from this Bayesian approach given 
the spectroscopic and photometric data from Johns-Krull et 
al. (f2008 ), the stellar evolution models from Girardi et al. (2002), 
and the orbit parameters determined from the Keplerian fit of 
the SOPHIE radial velocity measurements (§ 4. 1 1. The probabil- 



ity distribution function for the radius of XO-3 is centered near 
R* - 1.5 R , but extends with non-negligible density from 1.3 
to 2.0 R . It is well described by R+ = 1.6 ± 0.2 Rq. The corre- 
sponding masses are M* = 1.3 + 0.2 M G . This is a quantifica- 
tion of our "best guess" from the present observational data and 
prior knowledge about field stars. These parameters are reported 
in Table [2] 
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4.3. Transverse Rossiter-McLaughlin effect? 

The radial velocities of XO-3 measured with SOPHIE during the 
transit of January 28th, 2008 are plotted in Fig. [5] Surprisingly, 
they do not show the ordinary anomaly seen in case of prograde 
transit, i.e. a red-shifted radial velocity in the first half of the 
transit, then blue-shifted in its second half. During the full transit 
of XO-3b, the radial velocity is blue-shifted from the Keplerian 
curve, by about 100 ms _1 . Such shape is expected for a trans- 
verse Rossiter-McLaughlin effect, i.e. when the A misalignment 
angle is near 90° so the the planet crosses the stellar disk nearly 
perpendicularly to the equator of the star. This is apparently the 
case for XO-3b, whose transit seems to only hide some red- 
shifted velocity components, i.e. a part of the star rotating away 
from the observer. A schematic view of the XO-3 system with a 
transverse transit is shown in Fig. [6] 

We overplot in Fig. [5] models of Rossiter-McLaughlin ef- 
fects for XO-3b, for A = 0° (u pper p anel) and A = 90° (lowe r 
panel). Following Loeillet et al. d2008"l) and B ouchy et al. (2008), 
we used the analytical Ohta et al. (2005) description of the 
Rossiter-McLaughlin anomaly. We adopted the orbital param- 
eters of Table [2] a projected stellar rotation velocity Vsinl of 
18.5 kms -1 , and a linear limb-darkening coefficient e = 0.69 
from Claret d2004l >, for r eff = 6250 K and logg = 4.0 dex. 
The transit was centered on the t c time determined above from 
the February 29th, 2008 photometric transit. To take into ac- 
count for the large uncertainty in the masses and radii of the 
star and its planet derived from spectroscopic and light curve 
analyses (§[T]and § |4.2| i, we plot the models using two extreme 
sets of parameters over the SOPHIE radial velocities in each 
panel of Fig. [5} the solid line is the Rossiter-McLaughlin model 
with large masses and radii as favored from spectroscopic analy- 
ses, and the dashed line is the Rossiter-McLaughlin model with 
smaller masses and radii as favored by the light curve analysis. 
The dotted line is, for comparison, the Keplerian curve without 
Rossiter-McLaughlin effect. 

Table[3]summaries the parameters used for the different mod- 
els and the quantitative estimations of the quality of the fits. Note 
that the inclination used by Johns-Rrull et al. (2008) for large 
masses and radii, namely i = 79.32°, produces slightly too long 
transits duration when used together with our refined orbit. We 
used i = 78.6° in that case, which remains within the ±1.36° er- 
ror bar obtained on i by Johns-Rrull et al. (2008). Models with 
A = 0°, or without Rossiter-McLaughlin effect detection, pro- 
duce poor fits, with high^- 2 values and radial velocity dispersions 
of 60 to 75 ms _1 . This is significantly higher than the expected 
uncertainties on radial velocity measurements, around 33 ms~' 
( see T able [T| and the residuals of the Keplerian fit presented in 



Table|3j both produce similar fits with A = 90°. The;f 2 is slightly 
better in the case of large masses and radii but this does not seem 
to be significant according the noise level. 



§ |4.1| a(0 - C) = 29 ms" 1 . Thus, our SOPHIE data seem to 
exclude such ordinary solutions. 

The models with A = 90° produces lower with velocity 
dispersions of 42 or 44 ms -1 . The lower panel of Fig. [5] shows 
that transverse transits produce better fits of the data, centered 
on the expected mid-transit and with the adequate duration and 
depth. One should also note that the SOPHIE measurements per- 
formed just after the transit (orbital phases from 0.21 to 0.23) 
are well described by the Keplerian orbit model (see Figs. [2] 
and|5]l. We recall that these points were not used to determine the 



Keplerian orbit (§4.1 1; the orbital parameters were determined 
using only the first four measurements of the January 28th, 2008 
night (filled squares in Fig. |5), together with the measurements 
secured on other nights. The good match of the data with the 
A = 90° models argues for a transverse transit. This possible de- 
tection is independent of the set of stellar parameters adopted in 



E 



-1 1.4 



-1 1.6 



-12.0 



■1 1.6 



-1 1.8 



-12.0 







A = 0° J 


- 
- 


H 1 1 1 , 1 1 , 


^>>+ - 
— , — , — i — , — i ~> 


— ' — 1 — 




\ = 90° J 






tit, ■ 









0.16 



0.18 0.20 
phase 



0.22 



Fig. 5. Rossiter-McLaughlin effect models. Top: A = 0° (spin- 
orbit alignment). Bottom: A = 90° (transverse transit). On both 
panels, the squares (open and filled) are the SOPHIE radial- 
velocity measurements of XO-3 with 1-cr error bars as a func- 
tion of the orbital phase. Only the first four measurements (filled 
squares) are used for the Keplerian fit (together with 19 mea- 
surements at other orbital phases; see § 4.1 1. The dotted line is 
the Keplerian fit without Rossiter-McLaughlin effect. The two 
other lines show Rossiter-McLaughlin models with i - 78.6° 
and a/Ri, = 4.8 (solid line) and i = 84.9° and a/R^ = 7.2 
(dashed line). The summary of these parameters is in Table[3] 



Table 3. Chosen parameters for the Rossiter-McLaughlin effect 
models plotted in Fig.|5](see text). 
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The ~ 40 ms~' dispersion of the data from these transverse 
models remains slightly above the computed uncertainties on ra- 
dial velocity measurements. This suggests that some extra uncer- 
tainties might be present and not taken into account in the error 
budget. This make us considering this observation as a hint of 
detection for a spin-orbit misalignment. 

An explanation for this too large dispersion could be the 
high atmospheric refraction. Indeed, as seen in Sect. [2] the end 
of the transit was observed at large airmass. This could intro- 
duce biases in the radial velocity measurements that are difficult 
to quantify. This agrees with the increasing dispersion of the data 
from these transverse models, which is in the range 30-35 ms -1 
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Fig. 6. Schematic view of the XO-3 system with transverse tran- 
sit, as seen from the Earth. The stellar spin axis is shown, as 
well as the planet orbit and the A misalignment angle. The scale 
is in stellar radii. The limit between soft and strong grey on 
the /l-scale represents the favored value from our observations 
(A = 70°, see Sect. [5). 



in the first half of the transit, then in the range 40-45 m s _I in 
the second half. 

The larger dispersion might also be partly explained as the 
expected errors are increasing in the second part of the transit 
because of moonlight correction (see §[3}. In addition, it is pos- 
sible that the planet was crossing the stellar disk above a spot; 
this could cause extra radial velocity variations (jitter), as the 
anomaly that is visible near the phase 0.19. Stellar H and K 
Call lines do not show core emissions, but they are less deep 
than other F5 stars. This implies logR' HK = -4.6 + 0.2, and we 
can not exclude XO-3 presents stellar activity, including spots. 



5. A small radius for XO-3b 

Shortly after the submission of this paper, photometry of 13 tran- 
sits of XO-3b were released by Winn et al. (2008a). These new 
observations strongly favors the smaller values for XO-3 and 
XO-3b radii and masses. The parameters reported by Winn et 
al. d2008al) agree with the ones presented here (Table |2j; this 
strengthen the Bayesian approach we used in jj |4.2| Timing pa- 
rameters (as P or To) from Winn et al. (2008a) are more ac- 
curate thanks to their high-quality transit photometry, whereas 
orbital parameters (as e, to or K) are more accurate in the 
present study due to the high-quality radial velocity measure- 
ments with SOPHIE. 

The small R p value excludes a grazing transit for XO- 
3b, and the corresponding models plotted in solid lines in 
Fig. |5} the Rossiter-McLaughlin anomaly should thus be large 
and detectable, with an amplitude near the order of magnitude 
(Vsin/)(# p /fl*) 2 150ms" 1 (Winn et al. |2008a) . We fitted 
the SOPHIE data using the updated parameters from Winn et 
al. d2008ab . in particular R p = 1.217 R Jup , a/R+ = 7.07 and 
i = 84.20°. According^ 2 variations, the range of A compatible 
with our observations is 70° + 15°. The lowest^- 2 is ~ 64, im- 



plying a 42 ms~' velocity dispersions around the model. The 
best fit with these parameters is plotted in Fig. [7] The residuals 
are plotted in Fig. [8] in three cases: without transits, with spin- 
orbit alignment, and with A = 70°. Among them, the last case is 
clearly favored by our data when the parameters from Winn et 
al. J2008al ) are adopted. 
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Fig. 7. Rossiter-McLaughlin effect models with A — 70° and 
the small R p value reported by Winn et al. (2008a). The squares 
(open and filled) are the SOPHIE radial-velocity measurements 
of XO-3 with 1 -crerror bars as a function of the orbital phase. 
Only the first four measurements (filled squares) are used for 
the Keplerian fit (together with 19 measurements at other orbital 
The dotted line is the Keplerian fit without 



phases; see § 4.1 



Rossiter-McLaughlin effcet. The solid and dotted lines are the 
models with and without Rossiter-McLaughlin effect. 
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Fig. 8. Residuals of the Rossiter-McLaughlin effect fits. Top: 
Without transit. Middle: A = 0° (spin-orbit alignment). Bottom: 
A = 70°. The squares (open and filled) are the SOPHIE radial- 
velocity measurements of XO-3 with 1 -crerror bars as a func- 
tion of the orbital phase. Only the first four measurements (filled 
squares) are used for the Keplerian fit (together with 19 measure- 
ments at other orbital phases; see §4.1). The vertical, dashed line 
shows the center of the transit. 
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6. Conclusion and discussion 

Table [2] summarizes the star, planet and orbit parameters of the 
XO-3 system that we obtained from our analyses. The radial ve- 
locity measurements that we performed with SOPHIE during a 
planetary transit suggest that the spin axis of the star XO-3 could 
be nearly perpendicular to the orbital angular momentum of its 
pla net XO -3b (A = 70° ± 15°). We note that one Johns-Krull et 
al. (2008) HET measurement was obtained near a mid-transit of 
XO-3b. This radial velocity is blue-shifted by (260 + 194) m s" 1 
from the Keplerian curve, in agreement with the possible trans- 
verse Rossiter-McLaughlin effect we report here, though with a 
modest significance. 

The SOPHIE observation remains noisy, showing more dis- 
persion around the fit during the transit than at other phases. We 
consider this result as a tentative detection of transverse transit 
rather than a firm detection. Indeed, the end of the transit was ob- 
served at large airmasses, which could possibly biases the radial 
velocity measurements. Our fits favor a transverse transit, but 
one can not totally exclude a systematic error that would mimic 
by chance the shape of a transverse transit. This would imply that 
the radial velocities measured during the end of the transit night, 
at large airmasses, would be off by about 100 ms _1 , i.e. three 
to four times the expected errors. Other spectroscopic transits of 
XO-3b should thus be observed. They will allow the transverse 
Rossiter-McLaughlin effect to be confirmed or not, and to better 
quantify its parameters, such as the value of the misalignment 
angle A. 

Narita et al. (2008 ) estimate that the timescale for spin-orbit 
alignment through tidal dissipation is longer than thousand Gyrs. 
This timescale is uncertain, but much longer than the timescale 
for orbit circularization, which itself is longer than the age of 
XO-3, estimated in the range 2.4 - 3.1 Gyrs (Johns-Krull et 
al. 1200 8 ): there are thus no obvious reasons to exclude an eccen- 
tric, transverse system. A strong spin-orbit misalignment would 
favor of formation scenarii that invokes planet-planet scattering 
(Ford & Rasio 2006 ) or planet-star interaction in a binary system 
(Takeda et al. 2008 ) rather than inward migration due to inter- 
action with the accretion disk. This suggests in turn that some 
close-in planets might result from gravitational interaction be- 
tween planets and/or stars. Chatterjee et al. (2007 ) and Nagasawa 
et al. (120081 1 have recently shown that scattering with at least 
three large planets can account for hot Jupiters and predicts high 
spin-orbit inclinations (see also Malmberg et al. 2007 ). On an- 
other hand, XO-3b is an object close to the higher end of plan- 
etary masses. As discussed for instance by Ribas & Miralda- 
Escude (2007), there are some indications that these objects 
be low-mass brown dwarfs, formed by gas cloud fragmentation 
rather than core accretion; so that XO-3b may not necessarily 
constrain planet formation scenario. 

Finally, pseudo-synchronization might be questioned in the 
case of the massive XO-3b (Mxo-3 — 100 x Mxo-3bX which 
moves on an eccentric orbit with a periastron particularly near 
its host star. Tidal frictions might be high enough to tune the 
stellar rotation velocity close to the velocity of its companion 
on its orbit at the periastron (Zahn lT977b . The expected pseudo- 
synchronized stellar rotation is given by V lot = V p x a( f* e) , where 

V p = 2iij Jjz^ is the planet velocity at the periastron. For the 

XO-3 system, this translates into V rot =s SOg^-kms -1 accord- 
ing to the values in Table [2] As the XO-3 radius is larger than 
1.1 Rq, its rotation velocity Vsinl = 18.5 kms -1 is clearly 
smaller than the pseudo-synchronized velocity. However, we 
note that a spin-orbit misalignment would tend to reduce the 



pseudo-synchronized rotation velocity of the star. In that case, 
the planet approaches at the nearest of its star at low stellar lat- 
itude, and not above the stellar equator. Pseudo-synchronization 
might thus be possible if actually there is a significant spin-orbit 
misalignment in the XO-3 system. 
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